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Abstract. Christ and Kiselev J^,|| have established that the generalized eigenfunc- 
, tions of one-dimensional Dirac operators with L p potential F are bounded for almost all 

energies for p < 2. Roughly speaking, the proof involved writing these eigenfunctions as 
a multilinear series J2 n T n (F, ... ,F) and carefully bounding each term T n (F, ... ,F). It 
is conjectured that the results in || also hold for L 2 potentials F. However in this note 
^ | we show that the bilinear term T2(F,F) and the trilinear term Ts(F,F,F) are badly 

behaved on L 2 , which seems to indicate that multilinear expansions are not the right 
en ! tool for tackling this endpoint case. 

!: 

I. Introduction 

Let F(x) be a real potential on R. For each energy k 2 > we can consider the Dirac 
generalized eigenfunction equation 



in 



o 



on R. This Dirac equation can be thought of as a Schrodinger equation with potential 
V = F' + F 2 . For each k there are two linearly independent eigenfunctions <fi — (fik- A 
natural question from spectral theory is to ask whether these eigenfunctions are bounded 
(i.e. are in L^°) for almost every real k. In Christ and KiselevQ showed among other 
things that this was true when F 6 for any 1 < p < 2. It is well known (see e.g. 



12]) that the statement fails when p > 2, but the p = 2 case remains open. In M it is 



> 

^ ■ shown that for L 2 potentials one has absolutely continuous spectrum on [0, oo), but this 

is a slightly weaker statement. 

We briefly outline the arguments in @,||. The method of variation of constants sug- 
gests the ansatz 

<j){x) = a(x)e ikx + b(x)e~ ikx 

(-4- + FU(x) = -ika(x)e ikx + ikb(x)e- ikx . 
ax 

Substituting this into the previous and simplifying, we reduce to the first-order system 

a'{x) = F{x)e~ 2ikx b{x) 
b'{x) = F{x)e 2ikx a{x). 

For simplicity we may assume F is supported on the positive half axis. If we set initial 
conditions a(— oo) = 1, b(— oo) = for instance, and then solve this system by iteration, 



^^The results cited are phrased for Schrodinger operators but also extend to the slightly simpler case 
of Dirac operators, see B . 
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we thus obtain the formal multilinear expansions 

o = l+ T n (F,...,F); b= T n (F,...,F) 

"> 2 - even n>l, Odd 

where for each n > 1, T n is the n-linear operator 

T n (F h ... ,F n )(k,x) := I e~ 2ik ^U(- i y^F 1 (x 1 )...F n (x n ) d Xl ...dx n . 

J Xl<...<X n <X 

For integrable Fj we can define the n-linear operators 

T n (F u ... ,F n )(k, +00) := ! e- 2ik ^7=i(-^F 1 (x 1 )...F n (x n ) dx 1 ...dx n . 

J X\<...<X n 

The strategy of Christ and Kiselev was then to control each individual expression T n 
on L p . Specifally, they showed the estimate 

|| sup |T n (F, . . . ,F)(k,x)\\\ < C Pin \\F\\l Px (1) 

X k 

for all n > 1 and 1 < p < 2, where C p .„ was a constant which decayed rapidly in n and 
1/p + 1/p' := 1. In particular one has the non-maximal variant 

||T n (F,... +oo)|| . pV „,oc <C Pjn \\F\\ n LP . (2) 

The boundedness of eigenfunctions for almost every k then follows by summing these 
bounds carefully. 

It is tempting to try this approach for the endpoint p = 2. For n = 1 we see that 
Ti(F)(k, +00) is essentially the Fourier co-efficient F(k), while sup x \Ti(F)(k,x)\ is essen- 
tially the Carleson maximal operator CF(k). The estimates (0), ([!]) for p = 2 then follow 
from Plancherel's theorem and the Carleson-Hunt theorem respectively. 

For n = 2 the expression T 2 (F, F)(k, +00) is essentially H_(\F\ 2 )(k), where H_ is the 
Riesz projection 

H-F := X(-oo,o]F, 

and so (Q) follows for p = 2 by Holder's inequality and the weak- type (1,1) of the Riesz 
projections. We also remark that if the phase function x\ — x 2 in the definition of T 2 
were replaced by a\X\ +a 2 x 2 for generic numbers a±, a 2 then the operator is essentially a 
bilinear Hilbert transform and one still has boundedness from the results in 0, ||, |L3|. 

It may thus appear encouraging to try to estimate the higher order multilinear operators 
for L 2 potentials F. However, in this note we show 

Theorem 1.1. When p = 2 and n = 2, the estimate (|1|) fails. When p = 2 and n = 3, 
the estimate (0) fails. 

Because of this, we believe that it is not possible to prove the almost everywhere 
boundedness of eigenfunctions for Dirac or Schrodinger operators with L 2 potential purely 
by multilinear expansions; we discuss this further in the remarks section. 

The counterexample has a logarithmic divergence, and essentially relies on the fact 
that while convolution with the Hilbert kernel p.v .- is bounded, convolution with 

r X X 

or X( — ol(x) is not. It may be viewed as an assertion that L 2 potentials create significant 
long-range interaction effects which are not present for more rapidly decaying potentials. 
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Interestingly, our counterexamples rely strongly on a certain degeneracy in the phase 
function ^-Y x j on the boundary of the simplex x\ < . . . < x n . If one replaced this 

phase by Y^jXj, then we have shown in |§, 13 that the bound (§) in fact holds when 
p = 2 and n = 3. Indeed this statement is true for generic phases of the form . ctjXj. A 
similar statement holds for ([I]) when p = 2 and n = 2 and will appear elsewhere. 

The first author was supported by NSF grant DMS 0100796. The second author is 
a Clay Prize Fellow and is supported by a grant from the Packard Foundation. The 
third author was supported by a Sloan Fellowship and by NSF grants DMS 9985572 and 
DMS 9970469. The authors are grateful to M. Christ for pointing out the importance of 
the degeneracy in the phase function ^Y x j anc ^ f° r suggesting numerous valuable 

improvements to the manuscript. 



2. Proof of Theorem |1.1| 



The letter C may denote different large constants in the sequel. To be consistent with 
the previous notation we shall define the Fourier transform as 



We let iV ^> 1 be a large integer parameter, which we shall take to be a square number, 
and test (jl]), (fj) with the real- valued potential 

2N 

F(x):=J2 F i(x) 

j=N 

where the Fj are given by 

Fj(x) :=Ar-i C os(2^x)0(^-j), 

is a smooth real valued function supported in [—7,7] with total mass f (p = 1 such 

that 4> stays away from in [—1, 1], and A is a sufficiently large absolute constant whose 
purpose is to ensure that 

4 \^-Aj)\<m)\ 

jeZ\{0} 

for £ G [—1,1]. Informally, F is a "chirp" which is localized in phase space to the region 

{(k,x) : k = + O(i); x = Nj + O(N), N<j< 2N}. 

We may compute the Fourier transform of the Fj using the rapid decay of as 

Fj{k) = ^ e - 2i(Nk - Af) ^(Nk - Aj) + 0(N- 200 ) (3) 

in the region ^ < k < 3 A. We remark that the error term 0(iV~ 200 ) has a gradient which 
is also 0(iV- 200 ). 

Clearly we have ||i^||2 = 0(iV -1 / 2 ), and hence that 

11*112= 0(1). 
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We now compute 



in the region 



T 2 (F,F)(k,x) = I e 2ik{xi - X2) F{x 1 )F{x 2 ) d Xl dx 2 (4) 

' Xl<X2<X 

\Nk-A Jo \<l; x = N{j q -Vn+ 1 -) (5) 

for some integer ^ < j < 2N. In this region we show that 

\T 2 (F,F)(k,x)\> C' 1 log N, (6) 

which will imply that 

\\mp\T 2 (F,F)\(k,x)\\ L a,<*, > C-HogN 

x k 

and thus contradict (|]) for n = 2 and p = 2 by letting N go to infinity. 

We now prove (|6|). Fix k, j , x. Observe from (|j) that T 2 (Fj, Fj/)(k,x) vanishes unless 
j < j' < jo — y/N. Thus we may expand 

T 2 (F,F)(k,x)= HF^F^x) (7) 

N<j<j -VN 

+ T{F v F r ){k lX ). (8) 

N<j<j'<jQ-\fN 



We first dispose of the error term (8). In the region j < j 1 < j — viV, the conditions 
x\ < x 2 < x in (H) become superfluous, so we may factor 



T 2 (F j ,F f )(k,x) = F j (k)F j ,(k). 

However, since is rapidly decreasing and \j — jo\, \f — jo I > VN, we see from (H) that 

1^)1,14(^)1 <CN~ W0 . 

Summing this, we see that the total contribution of (8) is 0(iV~ 198 ). 
Now we consider the contribution of (7). We use the identity 

T 2 {F 3) F 3 ){k lX ) =T 2 {F 3) F 3 ){k 1+ ^) =H4\F J \ 2 )(k) (9) 

combined with The operator is a non-trivial linear combination of the identity 
and the Hilbert transform, while \F 3 \ 2 is essentially a non-negative bump function rapidly 
decreasing away from the interval [jA/N — 0(l/N),jA/N + 0(l/N)]. Because of this we 
see that for j ^ j we have 

H4\F 3 \ 2 )(k) = + 0(\j - j o n (10) 
J ~ Jo 



where c is a non-zero absolute constant. Summing this over all j < jo — yN and observing 
that j — jo has a consistent sign we see that the contribution of (7) has magnitude at 
least C- 1 \ogN, and (§) follows. 

We now compute T 3 (F, F, F)(k, +oo) in the region 

\Nk - Aj \ < 1; 1.4iV < j < 1.6JV. (11) 
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We will show that 

\T 3 (F,F,F)(k,+oo)\ ^C-MogiV (12) 

in this region, which will disprove (0) for n = 3 and p = 2 similarly to before. 

It remains to prove (|T^). Fix j . Observe that T 3 (Fj, Fj>,Fjti) (k,+oo) vanishes unless 
j < j' < j" . Thus we can split 

T 3 (F,F,F)(k,+oc) = T,{F v F v F j ){k,+^) (13) 

N<j<2N 

+ Y, T z {F v F J) F r ){k,+^) (14) 

N<j<j'<2N 

+ Yl r 8 (^ l F ij F J -)(*,+oo) (15) 

N<j'<j<2N 

+ T 3 (Fj, Ff, Fj")(k, +oo). (16) 

A r <J<i'<j"<2A r 

We first consider (13). We expand 

T 3 (Fj, F^F^ik, +oo) = j e^^+^Fjix^Fjix^F^xs) dx x dx 2 dx 3 . 

This is a linear combination of eight terms of the form 

N~ 3 / e 2ik{xr-x^) e 2i^{± X1 ± X2 ± Xi )^^l _ -^^2 _ _ ^ dxidx2dx3 . 

J Xi<X2<X% N N N 

making the substitutions y s := % — j for s = 1, 2, 3, this becomes 



'2/l<2/2<2/3 

for some phase e ld depending on all the above variables. 

We shall only consider the choice of signs (— y\ + y 2 — y 3 ); the reader may easily verify 
that the other choices of signs are much smaller thanks to stationary phase. In this case 
we can write the above as 

e ie I e 2i ^- A ^-^V(yi)0(y 2 )0(y3) d yi dy 2 dy 3 . 

Jy\<y2<y-i 

If kN — Aj = 0(1) we estimate this crudely by 0(1). Otherwise we can perform the y\ 
integral using stationary phase to obtain 



e" n ..,,J [ e^ kN - A ^<P(y 2 )<P(y 2 )<P(y 3 ) dy 2 dy 3 + 0(\kN - Aj]' 2 ) 



2i(kN - Aj) 



Performing another stationary phase we see that this quantity is 0(\kN — Aj\~ 2 ). Sum- 
ming over all j we see that (13) is 0(1). 

Let us now consider (16). When j < f < j", the constraints x\ < x 2 < x 3 in the 
definition of T 3 are redundant, and we can factorize 



T 3 (Fj, Fj,, F f ,)(k, +oo) = F j (k)F j ,(k)F jl ,(k). 



6 CAMIL MUSCALU, TERENCE TAO, AND CHRISTOPH THIELE 

Applying (|3|) and using the rapid decay of <p we see that 

^3(^,^,^)^+00)1 < C(l + \j-j \ + \j'-j \ + \j"-j \r 10 + CN- 100 . 

Summing over all we see that (16) is 0(1). 

It remains to control (15) + (14). First we consider (14). For this term the condition 
x 2 < x 3 is redundant, so we can factorize 



T 3 (Fj, Fj, F f )(k, +00) = T 2 (F 3 , Fj)(k, +oo)F r (k). 
Now consider (15). For this term the condition x\ < x 2 is redundant, so we can factorize 

T 3 (F f , F h Fj){k, +00) = F r (k) I e^-^F^Fjixs) dx 3 dx 2 . 



Writing x\ instead of x 3 we thus have 

T 3 (F f , Fj, Fj)(k, +00) = ^)(|F,(A;)| 2 - T 2 (F V F 3 )(k, +00)). 
Combining this with the previous, we thus see that 



(15) + (14)= ^(f-j)T 2 (F J ,F J )(k,+^)F f (k)+ £ F f (k)\F,(k)\ 2 . 

N<j,j'<2N N<j'<j<2N 



Using (H) as in (16) we see the second term is 0(1), so to prove flT^) it will suffice to show 
that 



I s S n 0' - j) T 2(Fj, Fj)(k, +00) Fj>(k)\ > _1 logiV. (17) 

N<j,j'<2N,j^f 

We first consider the terms with j' = jo- We claim these terms are the dominant contri- 
bution. From (pi), fliOD we conclude 



J2 sgnOo - 3)T 2 {F 3 , F 3 ){k, +oo)F jo (k) 

N<j<2N,j^j 

sgn(jo - j 



- E r^+^ oW + 0(1) . (18) 

Here c is the same non-zero constant as in (0), and Fj (k) is bounded away from by 
choice of <p. Thus the first term is greater than C~ l log N , so it suffices indeed to show 
that this term is the dominant contribution to (|i7|) . 

We consider the terms with j = j . Using that \T 2 (Fj, Fj)(k, 00) | < C we obtain 

£ \T 2 (F m ,F JO )(k,+oo)F~ik)\<C 

N<j'<2N,j ^j' 

This term is therefore negligible. 

Finally, we have to consider the terms with j, j' ^ jo- We have by the choice of A, 



£ \T2(Fj, Fj)(k, +oo)\\Fjt(k)\ 

N<j,j'<2N,j,j'^j 

A £ TJ+rTl^WI + C 

N<j<2N,j^j IJ JU| 
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This term is dominated by ([18]) . This completes the proof of ([12]) . 



3. Remarks 

• The counterexample can easily be extended to larger n (e.g. by appending some 
bump functions to the left or right of F). 

• The counterexample above involved a potential F which was bounded in L 2 , but for 
which sup.,, \T 2 (F, F)(k, x)\ and \T 3 (F, F, F)(k, +00) | were large (about logiV) on a 
large subset of [A, 2 A] . By letting N vary and taking suitable linear combinations 
of such variants of the above counterexample, one can in fact generate a potential F 
bounded in L 2 for which sup x \T 2 (F, F)(k, x)\ is infinite and \T 3 (F, F, F)(k, x)\ accu- 
mulates at 00 for x — > 00 for all k in a set of positive measure (one can even achieve 
blow-up almost everywhere). Thus it is not possible to estimate these multilinear 
expansions in any reasonable norm if one only assumes the potential to be in L 2 . 
Similarly if F had a derivative in L 2 ; it is the decay of F which is relevant here, not 
the regularity. 

• The unboundedness of T 3 on L 2 can be interpreted as stating that the (non-linear) 
scattering map F ^ bk{+oo) from potentials to reflection coefficients is not C 3 on 
the domain of L 2 potentials. Similarly the map F > a^+oo) from potentials to 
transmission coefficients is not C 4 on the domain of L 2 potentials. In particular 
these scattering maps are not analytic. 

• Despite the bad behavior of the individual terms Tk(F, ... , F), the transmission and 
reflection coefficients cik(x), bk(x) are still bounded for the counterexample given 
above. This phenomenon is similar to the observation that the function e lx = 1 + 
ix — x 2 /2 — . . . is bounded for arbitrarily large real x, even if the individual terms 
(ix) n /n\ are not. 

We now sketch the proof of boundedness of a^, b^. Suppose that k = Aj /N + 
0(1/N) for some N < j < 2N; we now fix j and k. We can write 

where G is the 2x2 matrix solving the ODE 

r»( \ ( F(x)e- 2ikx \ n( , n( , /l 0\ 

G (x) = I p^yikx q ) G{x); G(-oo) = I Q 1 J . 

We define the matrices Gj similarly by 

G'j(x) = ( Fj(x)e 2ikx 0^ ^ ) Gj(x); Gj{— 00) = ( ^ ^ ^ . 

We observe the identity 

G(x) = G jl {+oo)G jl - 1 (+oo) . ..G N (+oo) 

whenever N < j \ < 2N and x = N(ji + ^); this can be proven by an easy induction on 
ji and the observation that the above ODE are invariant under right-multiplication. 



s 
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One can compute the Gj (+00) using multilinear expansions (or using Gronwall's 
inequality), eventually obtaining 

/ 1 + J<L n \ 
G,(+oo)= ( >-* ^ +O(|i-i |- 2 ) 
V i-io / 

for all j 7^ jo, where C is a non-zero real constant. Because of the crucial factor of i 
in the diagonal entries we see that the operator norm ||Gj(+oo)|| of Gj is 



G,(+oo)|| = 1 + 0(|j-j'o 



-2<\ 



This allows one to multiply the Gj(+ 00) together and obtain boundedness of G(x) 
and hence cik(x), bk(x). 

In analogy with the observation concerning e tx , one may need to use the fact that 
F is real in order to obtain boundedness of eigenfunctions in the L 2 case. When F 
is real there are additional estimates available, such as the scattering identity 



log |a fc (+oo)| dk = C J \F(x)\ 2 dx 

for some absolute constant C; see for instance ||. 

We do not yet know how to obtain boundedness of eigenfunctions for L 2 potentials 
F. However we have been able to achieve this for a model problem in which the 



Fourier phases e are replaced by a dyadic Walsh variant e(k, x). See |]TT 



One can modify the counterexample to provide similar counterexamples for Schrodinger 
operators —-j-2 + V with V € L 2 , either by using the Miura transform V = F' + F 2 
mentioned in the introduction, or by inserting the standard WKB phase modification 
to the operators Tk as in 0. We omit the details. 

The multilinear expansion of a leads to an expansion of \a\ 2 , whose quadratic term 
is equal to 

2Re(T 2 (F,F)) = 2Re(H_(\F\ 2 )) = \F\ 2 

This term is in L 1 , which is better than the term T 2 (F, F), which is in general only in 
the Lorentz space L 1,oc . The higher order terms of the expansion of \a\ 2 are however 
unbounded again. Using the identity \a\ 2 = 1 + \b\ 2 we see that the fourth order 
term of \a\ 2 is equal to 

2Re(J\{F)T 3 (F, F, F)) 
We now define the modified potential 

G(x) = F(x)+G (x) 

where F is as in the proof of Theorem |1.1| and Gq(x) = <\>{x — A^ 3 ). Expanding the 
fourth order term by multilinearity, one observes that all terms can be estimated 
from above nicely with the exception of 



2Re(T 1 (G)T 3 (F,F,F)) 

Since T\[G) = G has more rapidly changing phase than Ts(F,F,F), the real part 
and the modulus Ti(G)T 3 (F, F, F) are of comparable size on a large set, and so this 
term is of the order log(iV) on a large set just like T 3 (F, F, F) itself. 
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